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Abstract 

We calculate the 0{g'^a) mixing coefficients of bilinear quark operators in lat- 
tice QCD using a standard perturbative evaluation of on-shell Green's func- 
tions. Our results for the plaquette gluon action are in agreement with those 
previously obtained with the Schrodinger functional method. The coefficients 
are also calculated for a class of improved gluon actions having six- link terms. 
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Typeset using REVT^ 



I. INTRODUCTION 



Symanzik's improvement program HTJ applied to on-shell quantities p| attempts to elim- 
inate cut-off dependence order by order by an expansion in powers of the lattice spacing a. 
To 0{a) in lattice QCD, this requires adding the 0{a) "clover" term to the Wilson quark 
action [Q. Quark operators also have to be modified by 0{a) counter terms, which generally 
involve new operators of higher dimension . In perturbation theory, the tree-level value 
of the clover coefficient and those of the counter terms of quark operators can be easily de- 
termined. They are sufficient to remove terms of 0{g^a\oga) in on-shell Green's functions 
evaluated at one-loop order, as has been explicitly demonstrated in Ref. To remove 
0{g'^a) terms which still remains, the counter term coefficients for quark operators have to 
be corrected by 0{g'^a) terms. For bilinear quark operators, these coefficients have been 
calculated in Refs. using the Schrodinger functional technique. 

In this article we analyze the 0{g'^a) coefficients of bilinear quark operators through 
a standard perturbative treatment of on-shell Green's functions of the operators. One- 
loop amplitudes with external quarks on the mass shell are expanded in powers of qa and 
rriRa, with q the external momenta and mn the renormalized quark mass, which leads to 
an alternative determination of the coefficients. Applying the procedure for the standard 
plaquette gluon action, we obtain results which are in agreement with those of Refs. JT],^ . 

Another application of our procedure is a calculation of the coefficients for gluon actions 
improved by an addition of six-link loop terms to the plaquette action. We treat three cases: 
the action |^ which is tree-level improved in Symanzik's sense to C(a^), and two types of 
actions |IO|JTT| improved by a renormalization-group treatment. The results should be useful 
in simulations employing improved actions for gluons. 

This paper is organized as follows. In Sec. II we define bilinear quark operators examined 
in this article, and write down renormalization relations between the renormalized and bare 
operators. Analysis of one-loop amplitudes are carried out in Sec. III. Numerical results for 
the coefficients and a comparison with previous work are made in Sec. IV. We close with 
some concluding remarks in Sec. V. 



II. CLOVER QUARK ACTION AND BILINEAR QUARK OPERATORS 

Consider the clover quark action defined by 

n 

- cswa^J2^^9ji'n(^^iuPf,A'n')'4'n- (II-l) 

We wish to construct a renormalized bilinear quark operator of form 

O^R = (^crV'c)^ , r = 1, 75, 7^, 7^75, a^^ (II.2) 

which is improved to 0{a) to one-loop order of perturbation theory, i.e., on-shell matrix 
elements of the operator have no errors of 0{a), 0{g'^a\oga) or 0{g'^a) when a — > with 
external momenta and the renormalized quark mass m^j kept fixed. 
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Our starting point is the tree-level improved operator of form, 



where the rotated quark fields ipc and -ip^ are given by 



ar f -* 
1 - - (i^^D^, 



ar 

1 

4 



mo 



mo 



with the covariant derivative defined as 

1 



U^{x)ijj{x + fi) — Uj^{x — jl)ijj{x — /t) 
ij{x + fi)Ul{x) - i){x - jl)U^{x - fi) 



(11.3) 

(11.4) 
(11.5) 

(11.6) 
(11.7) 



It has been demonstrated in Ref. that this operator is on-shell improved to 0{a) 
and 0{g^a\oga) with the tree-level value of the clover coefficient csw = 1. It has been 
noted furthermore that the field rotation can be generalized by using the equation of motion 
{p+ mo)^ = to 



ar 



z'^^,D^, - (1 - z)mQ 



^ ar 

~ y 



-zl^,D^, - (1 - z)mQ 



where 2; is a parameter. We then consider a generalized operator given by 

Ol = [l + ar (1 - z) mo] + zF® - z^V , 
where F*^ and F' are 0{a) and 0{a^) vertices defined as 



ar ( 



2 



F' 
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(11.8) 
(11.9) 

(11.10) 

fn.ii) 



The one-loop relation expected between the bare operator 
improved operator ( [11. 2|) has the form 

Ol = Z^^Ol - g^CpaniR Br - g^Cpa Cr O^, 



(11.12) 

11.101) and the renormalized 

(11.13) 



where Cp denotes the second-order Casimir eigenvalue for the quark filed, and the last 
two terms are needed to remove 0{g'^a) errors from on-shell matrix elements, with O^^ a 
dimension 4 operator with derivative. In a previous paper [0 we have evaluated for a 
class of improved gluon actions. Our task now is to generalize the analysis to (i) check that 
there are no 0{g^a\oga) errors with the operator ( [11.10| ), and (ii) determine the 0{g'^a) 
coefficients By and Cr- In the following we set the Wilson parameter r = 1. 



III. ANALYSIS OF ONE-LOOP AMPLITUDES 
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FIGURES 
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FIG. 1. Structure of one- loop diagrams for bilinear quark operator ( II.1C| ) 

The structure of one-loop diagrams relevant for our analysis is depicted in Fig. |I| where 
we indicate our momentum assignment. We calculate the corresponding amplitudes in Feyn- 
man gauge imposing the on-shell condition to external momenta, which becomes for these 
diagrams, 

(-z/ + mfi)r = r(^^ + mfi) = 0. (III.l) 

We identify the renormalized mass m/j, to be explicitly defined below, with the on-shell 
mass as there appear no divergences in these diagrams which should be renormalized into 
quark mass. We note that the bare quark mass mo enters in the field rotation ( [il.9| ) of the 
bare operator, where we make use of the tree-level equation of motion. 

The vertex function in momentum space calculated to one-loop order has the form 

= [1 + arriQ (1 - z)]T + za^ {ifT - Tif) 

+g^CF [1 + arrio (1 - z)] Tr + g'^CpzTr^ - g^Cpz'^Tr', (III.2) 

where F represents the tree level contribution, Tp, Tp® and Tp/ are one- loop contributions 
from the vertices F, F*^ and F'. 
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A problem in extracting 0{a) terms of the one-loop contributions is that they are infra- 
red divergent for on-shell external momenta. We treat this problem by supplying a mass A 
to the gluon propagator. The one-loop amphtudes, being functions of pa,p'a,miia and Aa, 
are then finite, which we expand around pa = p'a = and m^a = 0, keeping Aa finite. For 
this procedure to be justified, infra-red divergences which remain in the vertex function Gr 
after wave function renormalization should coincide with those in the continuum. We check 
this point explicitly below. 

As a first step to extract 0{a) terms, we expand one- loop contributions in terms of 
external momenta. Under the on-shell condition ( |111.1D the 0{a) term can be written in two 
alternative forms, i.e., 

T(r,r's>,r') = ^(r,r®,r')^ + '^(r,r®,r')^'^^P^^P^' (111.3) 

where g^^ = p'^ ± Pf^, and Q'^r^ is a dimension-4 operator vertex with the same quantum 
number as F as listed in Table | for each F. The on-shell identities which relate the two 
forms are given by 

-ic[^l5 + iqu(^iiul5 + 2mfi7^75 = 0, (III.5) 
-ia^uqt + i% - 2mi{7^ = 0, (III.6) 
?;7m75 = g^7/. = 0, (III. 7) 

+ 2imH75 = 0, (III. 8) 

g;7/. + 22mR = 0, (III.9) 
T^p^^ypiTiv - Opy-ip) - PpiripO^p - ^yOpp) + {ctplv - qtlii) + 4imRa^^ = 0. (III.IO) 

We choose to work with 5'^F+ and drop the + suffix. The momentum g+ = represents 
the momentum transfer at the operator vertex. We observe from the identities above and 
Table | that (J^F^ actually vanishes by the on-shell condition for the scalar and pseudo 
scalar operators. Substituting the expansion ( [1II.4|) into ( [1II.2| ), we obtain 



= [1 + aniQ {1 — z) + zamji] F + g'^Cp [1 + arrio (1 — z)] {VrT + v^iaqpT p) 

+g'^CFz{VY<s>T + vrrsiaq^f^) - g'^Cpz^iVYiT + vr'iaq^f^). (III. 11) 

As a second step, we expand the one-loop amplitudes V(r,r®,r') in terms of the lattice 
spacing a multiplied by the renormalized quark mass ttir. This leads to 

Vr = + \/r^°) + amnV^'^ + 0{a^), (III.12) 

Vr^ = h'^{T)amRL + V^""^ + amnV^^^ + 0{a^), (III.13) 
= V^r?) + amnV^}^ + 0{a^), (III.14) 



where V^p p/) are constants independent of a and g. In a similar expansion of f(r,r®,r')) 
we only need to keep the leading term in a, and hence they can be regarded as a constant 
as well. The logarithmic divergence L is defined as 
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with A being gluon mass to regularize the infra-red divergence which appears in on-shell 
vertex functions. The coefficients h2{T) and h2{T) are given in Table 0. 

In order to relate bare operators to renormalized operators we further need the renor- 
malization factor for quark wave function and mass Zm, which are defined as 

^0 = Z^'^'^n, (111.15) 

mo = Z-'mn + 9^CF—. (III.16) 

a 

The explicit form of these factors are obtained from the inverse full quark propagator ex- 
panded to 0{g^a), 

Sp^ = + mo + ^ap'^ -E{p,mo), (III.17) 

where the one-loop correction to self energy can be written as 

"So 



S(p,mo) = g Cf 



+ {-L + El) + mo (-4L + Eg) 



+ap'^ (-L + (Ti) + amo^^ (4L + (T2) + ami (2-^ + (^3] 



+0{a^) (III.18) 



with So,i,2 and cti 2,3 being constants. In terms of these constants, the renormalization factors 
are given by 

Z^^ = 1+ g^Cp (-L + El) + a m (-1 + g^Cp (L + 2 ai + (T2 - 3 Ei + E2)) , (III.19) 
= 1+ /Ci.(-3L-Ei + E2)+ am(^^ + ^^(-3L-2ai-2(T2 + 2(T3 + Ei)^ . (III.20) 

Here we have defined 

m = mo- g^Cp — , (III.21) 
a 

and have made use of the relation g^rrio = g^m + 0{g^). The wave function renormalization 
factor can be rewritten in terms of the renormalized mass, 

Z^^ = 1 + g^Cp i-L + El) + amR (-1 + g'Cp [aL - z„ + eS^^)) , (III.22) 

where 

Zm = -El + E2, (III.23) 

eS') = 2(Ti + (T2-3Ei + E2. (III.24) 

Replacing the quark mass tuq with the renormalized mass mn, we obtain for the the 
vertex, 
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1 + arriR + g^CpL 



l-z) + h'^iV)z - 3(1 - z) 



r 



+g^Cp [So(l -z) + V^'^ + zVi^'^ - z^vi^^ 
+g^CFamR - + (1 - z)V^°^ + V^'^^ + zV^ 



(1) 



1 + arriR + g^CpL 



'h2{T) 



r 



+g^CF\M^ ~ z) + V^'^ + zVl''^ 



r 



+^"CFamK [zm{l -z) + {l- z)V^°^ + V^'' + zV^ 
+(7^Cir(fr + zvr^ — z'^vr')iaq^T 



(1) 



(III.25) 



where the relation h2{T)/4: — /i2(r) = 3, vahd for each F, is used for the second equahty. 

We now muhiply the vertex function by the quark wave function renormahzation 
factor Z^-^ of ( [111.22| ). The 0{a) and 0{g'^aloga) terms all cancel out in the combination 
Z^^G^ for arbitrary values of the parameter z, and the result can be written as an operator 
identity ( [11.13| ) with the constants given by 

'h2{r) 



Zr' = l + g^CF 



1 ) L + Si + rr^°) + So(l -z)+ zV, 



(0) 



(III.26) 



Bt 



Si + S^') + V^'> - So(l -z) + z{-Zm - V^"" - V^^^^ + V, 



(0) 



Cr 



z'^vr' 



(III.27) 
(III.28) 



Let us add a remark on infra-red divergence. For f(r,r®.r') that enter in the coefficient 
Cr, these divergences cancel out in the total contribution at each order of z. On the other 
hand, infra-red divergences remain in V(r,r®,r')- The divergence in the 0{g'^a) term in G^ 
are, however, canceled by that of the quark wave function renormahzation factor. The 
remaining divergence, which is of form g^L and appears in Z-p, coincides with that which is 
present in on-shell vertex functions for the renormalized operator in the continuum. 

In the calculation above we employed the bare operator ([II.10|) which contains the bare 
quark mass mo. It is possible to replace mo by the subtracted mass m of ( |111.21| ). Defining 



= [1 + ar (1 - z) m] ^jT^j + zT^ 



z^T\ 



(III.29) 



it is straightforward to check that 0{a) and 0{g'^aloga) terms also cancel for this operator; 
So plays no role in improving the operator. The renormahzation coefficients for this operator, 
which is more convenient for practical use, are obtained from ( [III.26| - pI.28|) by eliminating 
So. 
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IV. RESULTS FOR THE ONE-LOOP COEFFICIENTS 



Manipulations in the previous section have reduced the determination of one-loop co- 
efficients to an evaluation of a number of integral constants. Working out the integrands 
for the integrals is a straightforward but tedious task, which we carry out by Mathematica. 
The output is converted to a FORTRAN code, and the integrals are evaluated by the Monte 
Carlo routine VEGAS in double precision. We employ 20 sets of 10^ points for integration 
except for Ca for z = for which we use 20 sets of 10^ points. Errors are estimated from 
variation of integrated values over the sets. 

For the gluon action we consider the form given by 

^gluon = I Co Trt/pZ + Ci ^T^Urtg + C2 Tff^chr + C3 Y Trf/p;^ 

3 y plaquette rectangle chair parallelogram 

(IV.l) 

where the ffist term represents the standard plaquette term, and the remaining terms are 
six-link loops formed by a 1 x 2 rectangle, a bent 1x2 rectangle (chair) and a 3-dimensional 
parallelogram. The coefficients Cq, • • ■ , C3 satisfy the normalization condition 

Co + 8ci + I6C2 + 8C3 = 1. (IV.2) 

At the one-loop level, the choice of the gluon action is specified by the pair of numbers 
Ci and C23 = C2 + C3. We calculate the constants for five cases: (i) the standard plaquette 
action ci = 0, C23 = 0, (ii) the tree- level improved action in the Symanzik approach ci = 
— 1/12,C23 = 0, and (iii) three choices suggested by an approximate renormalization- 
group analysis, ci = — 0.331, C23 = and Ci = —0.27, C23 = —0.04 by Iwasaki and 
ci = -0.252, C23 = -0.17 by Wilson fig. 

Let us write O^-^ for the renormalized bilinear operator in the continuum in the MS 
scheme, and define 

^175 = ^175^0, (IV.3) 
where Cq denotes the lattice bare operator and 

Results for the finite constant z-p for the lattice operator (p.l.lOD rotated with the bare mass 



mo have already been given in a previous paper . We list in Table |T| the values of z-p 



for the operator ( [1II.29|) defined with the subtracted mass m at z = for completeness. 

Our new results for the one-loop coefficients Cr are given in Table |V], and those for 
By in Tables and Vl for the definition excluding Eg. Numerical values are given for the 



coefficients of expansion in z defined as 

Br = BP+zB^^^-z'B^^\ (1V.5) 
Cr = C^'^ + zC^}^ - z^C^'l (1V.6) 
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In the tensor channel only Cj^^ and Bj^^ are evaluated as necessity for the operator in this 
channel does not seem to warrant a CPU time-consuming calculation of integrands which 
are more complex than the other cases. In Tables and ^ a general trend is apparent 
that the coefficients are reduced by roughly a factor two for renormalization-group improved 
gluon actions as compared to those for the plaquette action, as already observed for zr |]12[ . 

Comparison of our results with those of Refs. obtained with the Schrodinger func- 
tional is made in the following way. The authors of these references start from a local bilinear 
operator 



and relate it to the renormalized operator through 

0^,=Z^\g',a)Ol-Cr^q,^l. 
The renormalization factor is expanded in the lattice spacing a as follows, 

^r(/, a) = Zorig^) f 1 + ambr{g' 



(IV.7) 
(IV.8) 
(IV.9) 

where Zor(fi'^) does not contain terms of 0{a), 0{g'^aloga) or 0{g^a). The renormalization 
factor for the quark field ip and the quark mass m are expanded in a similar manner. 
With these definitions, the vertex function of the bare operator has the form 

= T + g^Cp (^^L + V^^^ + amV^^^^ T + g'^CFawiq^.f ^. (IV. 10) 

The expressions for the wave function, quark mass and quark bilinear operator renormaliza- 
tion factors are given by 



^1/2 



1 + 



g'c. 



1 

2 



1 + am ( — - + g'^Cp 



,0-2 V , ^2 



(l + /C^(3L + (Si-S2))) (l + 



am 



(Tl + 0-2 - 0-3 



Si + -S2 
2 



1 



g'Cp ( (M£) - 1 j L + Si + Vl'^^ ) (1 + am (1 - g'Cp (Si + S « + V^'^ 



Making an expansion 

bx 

cr = g'Cpc\." + ■ 
we find that the 0{g'^a) coefficients are given as 



5(0) +,2^^5(1) + . 
.(1) 



X 



T,^,m, 



5(1) 



, 0"2 V , ^2 

ai + y-Si + -, 



(Xi + 0-2 — (73 



Si + -S2m, 



Si + S« + 



-vr. 



(IV. 14) 
(IV.15) 



(IV.16) 

(IV. 17) 

(IV. 18) 
(IV.19) 



(IV.ll) 
(IV.12) 
(IV. 13) 
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for z 



Comparing these expressions with ([III. 27 ) and ( [III.2 
tabulated in Table and 6p ■* equals our Bf for z 



we see that Cp equals our C^^ 
without the Sq term given 



in Tables and Our results for 6^^^ and \))^ are given in table |VI1| , where we also list 
the contribution from the wave function renormalization factor = —(Si + T.[ 



-2b: 



(1) 



In Table [VIII| we collect our results for the mixing coefficients for z = for the plaquette 
gluon action, and compare them with those of Refs. [0,^. As we already remarked, we 
employ 20 sets of 10^ points for evaluating c^"* by VEGAS in this table, with which we find 
a complete agreement with the result of Refs. Good agreement is also found for all 

the other coefficients obtained with 20 sets of 10^ points. We do not pursue more precise 
evaluation for the latter coefficients since it would require significantly more computing 
power due to an increased complexity of integrands and the number of terms. 

It was observed in Ref. that the values of 6p ^ are close to each other. Numerically 
this arises from the fact that the contribution from the wave function renormalization Bq, 
common to various Dirac channels F, dominates over the vertex contributions. 



V. CONCLUDING REMARKS 

In this article we have carried out a perturbative evaluation of vertex functions to deter- 
mine the 0{g'^a) mixing coefficients of bilinear quark operators. For the standard plaquette 
action for gluons, our results agree with those obtained previously with the Schrodinger 
functional method. We have also generalized the determination to a class of improved gluon 
actions for use in numerical simulations employing such actions. 

Our calculations are carried out by an expansion of vertex functions regarding external 
momenta and renormalized quark mass m^j as small in units of lattice spacing a. Hence the 
present work does not cover the case of heavy quark such that m/ja > 0(1). It has been 
pointed out recently in a one-loop calculation in NRQCD ]16[ that the mixing coefficient 



c^'' for heavy-light axial vector current is large compared to the value for the light-light 



case treated here. In our calculation a significant cancellation is observed between terms 
from various diagrams contributing to c^"*. To understand whether the large value of c^'* 
for heavy quarks results from lifting of such a cancellation requires an extension of our 



calculation without making an expansion in m/ja 13-15 . 
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TABLES 



TABLE I. Dimonsicm 4 vertex mixing with ipTip at one loop order. 


r 


<r+ q-r- 


7a. 

75 
1 


5^7m75 Qillnlb 


TABLE IL Coefficients of logarithmically divergent term. 


r 


/i2(r) hm 


7m75 
7m 

75 

1 


4 -2 
4 -2 
16 1 
16 1 
-3 
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TABLE III. Finite part zr of renormalization factor for bilinear quark operators {z = 0). 
Coefficients of the term c'g^r{n = 0, 1, 2) in an expansion zr = Zy^ + csw^^^ + (^'sw^r^ given 
in the column labeled (n). Errors are at most in the last digit given. 



gauge action 

ci ('■2:i 


V 

((-)) (1) (2) 


(0) 


A 

(1) 


(2) 








-20.618 


4.745 


0.543 


-15.797 


-0.248 


2.251 


-1/12 





-16.603 


4.228 


0.464 


-12.540 


-0.198 


2.021 


-0.331 





-11.099 


3.326 


0.336 


-8.192 


-0.125 


1.610 


-0.27 


-0.04 


-11.540 


3.418 


0.353 


-8.523 


-0.131 


1.657 


-0.252 


-0.17 


-10.525 


3.248 


0.338 


-7.707 


-0.117 


1.587 



gang 


3 action 




S 






P 






T 




Cl 


C23 


(0) 


(1) 


(2) 


(0) 


(1) 


(2) 


(0) 


(1) 


(2) 








-12.953 


-7.738 


1.380 


-22.596 


2.249 


-2.036 


-17.018 


3.913 


1.972 


-1/12 





-9.607 


-6.836 


1.367 


-17.734 


2.015 


-1.745 


-13.539 


3.490 


1.719 


-0.331 





-4.858 


-5.301 


1.266 


-10.673 


1.601 


-1.281 


-8.939 


2.751 


1.300 


-0.27 


-0.04 


-5.260 


-5.454 


1.292 


-11.292 


1.644 


-1.316 


-9.283 


2.827 


1.344 


-0.252 


-0.17 


-4.366 


-5.166 


1.287 


-10.001 


1.565 


-1.212 


-8.427 


2.687 


1.271 



TABLE IV. Mixing coefficients Ca, Cy, Ct for axial vector, vector and tensor currents. 
Coefficients of the term z"'{n = 0, 1, 2) are given in the column marked as {n). C^^ and C^^ are 
not calculated. 



gauge action 

Cl C23 


Ca 

(0) (1) (2) 


Cv 

(0) (1) (2) 


Ct 

(0) 



-1/12 
-0.331 
-0.27 -0.04 

-0.252 -0.17 


-0.005680(2) -0.002316(3) -0.003808(1) 
-0.00451(1) -0.00179(2) -0.002781(7) 
-0.00285(1) -0.00099(2) -0.001579(6) 
-0.00302(1) -0.00108(2) -0.001660(6) 

-0.00281(1) -0.00098(2) -0.001131(6) 


-0.01226(3) -0.01560(4) -0.000217(9) 
-0.01030(4) -0.01262(4) 0.000034(8) 
-0.00729(4) -0.00825(4) 0.000175(7) 
-0.00757(4) -0.00858(4) 0.000159(7) 

-0.00705(4) -0.00772(1) 0.000153(7) 


-0.00898(1) 
-0.00741(1) 
-0.00508(1) 
-0.00530(1) 

-0.00195(1) 
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TABLE V. Mixing coefficients Ba, By, Bt for axial vector, vector and tensor currents. Co- 
efficients of the term z"'{n = 0, 1,2) are given in the column marked as (n). b!^^ and b!^^ are not 
calculated. 



gauge action 




Ba 






By 




Bj' 


Cl 


C23 


(0) 


(1) 


(2) 


(0) 


(1) 


(2) 


(0) 








0.1141(1) 


-0.0846(1) 


0.01637(3) 


0.1150(2) 


-0.0442(2) 


0.03255(5) 


0.1044(1) 


-1/12 





0.0881(1) 


-0.0666(1) 


0.01328(3) 


0.0886(2) 


-0.0353(2) 


0.02551(4) 


0.0795(1) 


-0.331 





0.0547(1) 


-0.0419(1) 


0.00867(3) 


0.0550(2) 


-0.0228(1) 


0.01583(4) 


0.0482(1) 


-0.27 


-0.04 


0.0572(1) 


-0.0438(1) 


0.00909(3) 


0.0575(2) 


-0.0238(1) 


0.01656(4) 


0.0505(1) 


-0.252 


-0.17 


0.0512(1) 


-0.0393(1) 


0.00827(3) 


0.0514(2) 


-0.0213(1) 


0.01479(4) 


0.0448(1) 



TABLE VL Mixing coefficients Bp, Bs for pseudo scalar and scalar density. Coefficients of 
the term z'^{n = 0, 1, 2) are given in the column marked as (n). 



gang 


e action 




Bp 






Bs 




Cl 


C23 


(0) 


(1) 


(2) 


(0) 


(1) 


(2) 








0.1148(1) 


-0.01698(6) 


0.03788(5) 


0.1444(2) 


-0.0535(2) 


0.00864(6) 


-1/12 





0.0890(1) 


-0.01378(6) 


0.02872(4) 


0.1144(2) 


-0.0403(2) 


0.00718(5) 


-0.331 





0.0561(1) 


-0.00909(6) 


0.01697(4) 


0.0747(2) 


-0.0235(1) 


0.00475(4) 


-0.27 


-0.04 


0.0586(1) 


-0.00937(6) 


0.01782(4) 


0.0777(2) 


-0.0248(2) 


0.00501(4) 


-0.252 


-0.17 


0.0527(5) 


-0.00816(6) 


0.01569(4) 


0.0706(2) 


-0.0221(2) 


0.00452(4) 



TABLE Vn. Mixing coefficients , bin for quark operator. Values for b^ are also included. 





^auge action 








Cl 


C23 














-0.05191(3) 


-0.07218(5) 


0.10381(9) 


-1/12 





-0.03968(3) 


-0.05722(5) 


0.07937(9) 


-0.331 





-0.02430(3) 


-0.03737(5) 


0.04860(9) 


-0.27 


-0.04 


-0.02543(3) 


-0.03891(5) 


0.05086(9) 


-0.252 


-0.17 


-0.02262(3) 


-0.03526(5) 


0.04525(9) 



TABLE VIIL Comparison with previous work for the plaquette gluon action. 





^A 


Cy 


6^'^ 


"a 




Op 


Ours 


-0.005680(2) 


-0.01226(3) 


-0.07218(5) 


0.1141(1) 


0.1150(2) 


0.1148(1) 


Sint-Weisz [| 


-0.005680(2) 


-0.01225(1) 


-0.07217(2) 


0.11414(4) 


0.11492(4) 


0.11484(2) 
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